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Abstract. In this paper we prove uniqueness for an inverse boundary value 
problem for the magnetic Schrodinger equation in a half space, with par- 
tial data. We prove that the curl of the magnetic potential A, when A G 
ompi^'^^^^)^ and the electric pontetial q E L^^^p(IR|^, C) are uniquely de- 
termined by the knowledge of the Dirichlet-to-Neumann map on parts of the 
boundary of the half space. 



1. Introduction 
We consider a magnetic Schrodinger operator La^q, defined by 

3 

(1.1) LA,g ■■^Y.i^ + ^j(^))' + 

on the half space M"?. := {a; S R'^ : < 0}. We assume that 

(1.2) AeW^^^piW,^'), ^ndqeLZnpiW,C), lmg<0. 
Here 

W^^piW,^^) ■■={A\^ ■■ Ae V^^l^°°(M^K^),supp(A) C compact} 
and similarly, we define 

L^^^p{W,C) {q e i°°(i^,C) : supp(9) C W. compact}. 
Consider the Dirichlet problem 

(La a - fc^)w = in R^ 
(1-3) I , 

where fc > is fixed and / e -ffcomp(9M^). Here we also require that the solution 
u should satisfy a boundary condition at infinity, which will be the Sommerfeld 
radiation condition 

(1.4) lim \x\(^^-iku{x)] =0. 

|x|-J-oo \ 0\x\ J 

Solutions satisfying this condition are called outgoing or radiating solutions. We 
will also occasionally use the term incoming solution. This refers to a solution of 
(|1.3p that satisfies (|1.4p . when the factor ~ik is replaced by ik. The existence and 
uniqueness of a solution u € Hf^^{M.^_) to the problem (|1.3p and (|1.4p is proven 
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in |19) . This allows us to define the so called Dirichlet to Neumann map A^^^, 
(DN-map for short), A^.g : H!/,lp{dM.l) Hf^^idRl) as 

/ ^ (dn + iA ■ ?i)w|aR3_, 

where u is the solution of the Dirichlet problem p. 3D . (II. 4[) and / is the value of u. 
Here n = (0,0, 1) is the unit outer normal to the boundary 9M'^. 

The inverse problem is to investigate if the DN-map uniquely determines the 
potentials A and q in R^. It turns out that the DN-map does not in general 
uniquely determine A. This is due to the gauge invariance of the DN-map, which 
was first noticed by [52]. 

Lemma 1.1. Let A e iyi'°°(R3, M^) ^ ^ L°°{W). Then 
(i) For all ip (E C1'1(M3,R) we have 

e-''f'AA,qe''f' = Aa+vv,9- 
(ii) There exists ip 6 C"'"'"'"(M'^, M) with 4'\{x3=o} — 0^ /o?^ which 

^A,q = AA+Vlf>,q 

and (A-h VV')|{:,3=o} (Ai,A2,0). 
Proof. See [ig. □ 

Part (ii) of this Lemma shows that A^.^ cannot uniquely determine A, since we 
can change a potential by a gauge transformation without changing the DN-map. 
The DN-map does however carry enough information to determine V x A, which 
is the magnetic field in the context of electrodynamics. 

When considering a pair of magnetic potentials Aj, j = 1, 2, we use the notation 
Aj = { Aj, I, Aj Aj ^3) for the component functions. Furthermore we let C 
be a relatively open and bounded set, for which 

y supp(Aj) U supp(gj) C A^, 

and for which d N piecewise and M.^_ \ N is connected. 

We now state the main result of this paper, which generalizes the corresponding 
results of [13] , obtained in the case of the Schrodinger operator without a magnetic 
potential. 

Theorem 1.2. Let Aj G W^i^^p(R3, r3) and qj G i^^p(R3,C) be such that 
Imqj <0,j~ 1,2. Let ri,r2 C (JK'l be open sets such that 

(1.5) (aM^ \]v)nr, ^0, J = 1,2. 

Then if 

(1-6) A^,,,,(/)|r, =A^,,,,(/)|r,, 

for all f e H^lpidRl), supp(/) C Fa, then 

V X Ai = V X A2 and qi = q2 in 

We would like to emphasize that in Theorem ll.2[ the set Fi, where measurements 
are performed, and the set F2, where the data is supported, can be taken arbitrarily 
small, provided that (jl.Sp holds. The result of Theorem 11.21 pertains therefore to 
the inverse problems with partial data. Such problems are important from the point 
of view of applications, since in practice, performing measurements on the entire 
boundary could be impossible, due to limitations in resources or obstructions from 
obstacles. 
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The first uniqueness result, in the context of inverse boundary value problems 
for the magnetic Schrodinger operator on a bounded domain, was obtained by Sun 
in [22], under a smallness condition on A. Nakamura, Sun and Uhlmann proved 
the uniqueness without any smallness condition in [17] . assuming that A S C^. 
Tolmasky extended this result to magnetic potentials in [53], and Panchenko 
to some less regular but small magnetic potentials in |18| . Salo proved uniqueness 
for Dini continuous magnetic potentials in j20j . The most recent result is given 
by Krupchyk and Uhlmann in fl2| . where uniqueness is proved for L°° magnetic 
potentials. In all of these works, the inverse boundary value problem with full data 
was considered. 

In [B], Eskin and Ralston obtained a uniqueness result for the closely related 
inverse scattering problem, assuming the exponential decay of the potentials. The 
partial data problem in the magnetic case was considered by Dos Santos Ferreira, 
Kenig, Sjostrand and Uhlmann in [5] and by Chung in [4j. 

The inverse problem for the half space geometry, without a magnetic potential 
was examined by Cheney and Isaacson in [5]. The uniqueness for this problem in 
the case of compactly supported electric potentials was proved by Lassas, Cheney 
and Uhlmann in |13| , assuming that the supports do not come close to the boundary 
of the half space. The result of Theorem 11.21 is therefore already a generalization 
of the work [13], even in the absence of magnetic potentials. Li and Uhlmann 
proved uniqueness for the closely related infinite slab geometry with ^ = 0, in [16] . 
Krupchyk, Lassas and Uhlmann did this for the magnetic case in [11] . In both of 
these works, the reflection argument of Isakov [5] played an important role. The 
uniqueness problem for the magnetic potentials in the slab and half space geometries 
has also been studied in a recent paper by Li [15]. The half space results in [TS] 
differ from the ones given in this work, by concerning the more general matrix valued 
Schrodinger equation and by assuming regularity on the magnetic potential. 

The half space is perhaps the simplest example of an unbounded region with 
an unbounded boundary. It is of special interest in many applications, such as 
geophysics, ocean acoustics, and optical tomography, since it provides a simple 
model for semi infinite geometries. We would like to mention that the magnetic 
Schrodinger equation is closely related to the diffusion approximation of the photon 
transport equation, used in optical tomography [T]. 

The paper is organized as follows. Section 2 contains a review of the construc- 
tion of complex geometric optics solutions for magnetic Schrodinger operators with 
Lipschitz continuous magnetic potentials. In section 3 we derive the central integral 
identity. The proof of Theorem 1 1.2 1 is contained in sections 4 and 5. The appendix 
contains an extension of Green's second formula and a statement of the unique 
continuation principle for easy reference. 

2. Complex geometric optics solutions 

Let 17 C be a bounded domain with C°°-boundary, and let A e W'^'°°{n, M^), 
q € L°°{il, C). The task of this subsection is to review the construction of complex 
geometric optics solutions for the magnetic Schrodinger equation, 

(2.1) LA,qU = in n. 

A complex geometric optics solution to (j2.ip is a solution of the form 

(2.2) u{x, C; h) = e-<'^{a{x, (; h) + r(z, C; h)), 

where C, € C'^, ■ = 0, a is a. smooth amplitude, r is a remainder, and ft. > is a 
small parameter. 
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In the case when A G (7^(51) and q £ L°°(ri), such solutions were constructed in 
[5] using the method of Carleman estimates, and the construction was extended to 
the case of less regular potentials in [9] , see also [11] . 

Let ip{x) = a ■ X, a € M.^, \a\ = 1. The fundamental role in the construction of 
complex geometric optics solutions is played by the following Carleman estimate, 

(2.3) hMnU^) < C\\e^/''h^LA,ge-'^/"u\\L2^n), 

valid for all u € C^{^1) and < h < ho, which was proved in [S| and [S]- Here 

\\u\\hI,{Q) = Wuh^Q) + \\hS7u\\L2^ny 

Based on the estimate (|2.3p . the following solvability result was established in 
P Proposition 4.3]. 

Proposition 2.1. Let A e W^'°°{n,R^), q G L°°(f7,C), a G M^, |a| = 1 and 

f{x) = a ■ X. Then there is C > and ho > such that for all h £ {0, ho], and any 
f G L'^{fl), the equation 

e'^'^'h^LA,qe-'^'^u = f m 
has a solution u G H^{VL) with 

C 

Our basic strategy in constructing solutions of the form (j2.2p is to write (j2.ip , 

as 

(2.4) L^r = -L^a, 
where ~ e-^^/'^/i^L^^e^ C/'i. rj.^^^ 

we first search for a suitable a, after which 
we will get r by Proposition 12.11 We must however take some care in choosing a 
and the way it depends on h, since we need later that Hr-H^i (q-j — > 0, sufficiently 
fast as ft, — > 0. We need a also to be smooth enough. This will be handled as in [9]. 

We extend A G VF^'°°(f2, M'^) to a Lipschitz vector field, compactly supported 
in f2, where f2 C is an open bounded set such that fl CC (l. Wc consider the 
mollification A" := G Cg°{0,,R^). Here e > is smaU and ipeix) = e'^-tp^x/e) 

is the usual moUifier with -ip G C^{W^), < ip < 1, and J -ijjdx = 1. We write 
A^ — A ~ AK Notice that we have the following estimates for A^ , 

(2.5) \\A'\U^^n)^Oie), 
I|a"A«|Uoo(o) =0(e-l"l) for all a, 

as e — >■ 0. 

We shall work with a complex C = Co + Ci depending slightly on h, for which 

(2.6) C • C = 0, Co := a + i^, a, /? G 5^ a • /3 = 0, 

Co independent of h and Ci = 0{h), as h ^ 0. 
By expanding the conjugated operator we write the right hand side of (|2.4p as 
L^a =^{-h^A - 2i{-iCo + hA) ■ /iV - 2Ci ■ hV + h^A^ 

(2.7) - 2ihCo ■ {A^ + A^) ~ 2ihCi ■ A - ih^{V ■ A) + h^q)a. 

Now we want a to be such that this expression decays more rapidly than 0{h), as 
h-^0. 

Consider the operator in (|2.7p . ignoring for the time being a and its possible 
dependence on h. We would like to eliminate from this operator the terms that are 
of first order in h. Notice first that Ci = 0{h) and that we can control ||^'' ||L=o(n) 
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with h, if we choose e to be dependent on h. Then in an attempt to ehminate first 
order terms in h, it is natural to search for an a for which 

(2.8) . Va = -<o ■ A^a, in fl. 

We will look for a solution of the form a = e*. The above equation becomes 
then 

(2.9) Co • V$ = -<o • A", infi. 

Pick a 7 6 5^, such that j-L{a, /?}. 

Next we consider the above equation in coordinates y, associated with the basis 
{a, j3, 7}. Let T be the coordinate transform y = Tx := (a; • a, x • /3, x • 7). Using 
the chain rule and the fact that T^^ = T* , one gets thalQ 

V($ o T-^){Tx) = r[V$(a;)]*. 

We therefore have that 

(1, i, 0) ■ V($ o T~^){Tx) = (1, i, 0) • r[V$(a;)]* 

= (a- V + i/3- V)$(x) 
= Co • V$(a;). 

Equation (|2.9p gives hence the 9-equation 

(2.10) 2 d, •($ o T-i)(y) = -zCo ■ (A« o r-i)(y), 

where dz — {dy^ -\-idy.^l2. We will solve this using the Cauchy operator 

N-^![x):^- f ^. fix~is,,S2,0))ds,ds2, 

which is an inverse for the 9-operator, N := {dy^ + idy^)/2 (see e.g. [7] Theorem 
1.2.2). We will need the following straightforward continuity result for the Cauchy 
operator. 

Lemma 2.2. Let r > and / G W''-°°{"M?), k > and assume that supp(/) C 
B(0,r). Then 

II^~"^/IIh"'-'°°(R3) < C'fe||/llw"=- = (R3) 

for some constant Cfe > 0. If / G Co(M.^), then N-^f e C(M^). 

Proof. See e.g. [21]. □ 

Returning to (PTU)) we get that $ = i7V-i(-<o ■ (^" oT^i)) oT. More explicitly 
we have 

<f,( , ^ 1 /■ -<o-^"(a^-J^"^(si,S2,0)) ^ , 

(2.11) <I>(.T,Co;e) = / dsids2, 

where T~^(si, S2, 0) = sia + S2/3. We have thus found a solution a = e* to equation 
(|2.8p . We will choose e so that it depends on h, which implies that a will depend 
on h. In order to determine how the norm of r will depend on h and also for 
later estimates, we will need to see how |j d" a||ioo depends on h. Lemma [2.21 and 
estimate p.5[) imply the following result. 



Lemma 2.3. Equation (j2.8l) has a solution a G C°°{Cl) satisfying the estimates 
(2.12) ||9"a|Uoo(o) <ae-l"l for all a. 



^Here T* is the transpose of T. 
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We can now write the L°°(ri) norm of (|2.7p as 

||ica||L°=(J2) = II - h'^LA.qa + 2ihCo ■ + 2Ci • hVa + 2ihC,i ■ Aa||ioo(Q). 

Using (P3|) . (pn^ and the fact that Ci = 0{h) we have that 

\\L^a\\L^(n) ^ 0{h^e~^ + he). 

Choosing e = h^/'^, gives finally ||L^a||x,=°(o) = 0{h'^/^), as /i ^ 0. 
Finally to solve (|2.4p for r, we rewrite it as 

(2.13) e-^-^'=«/''/i2i^,je^-R<=C//i(e«-imC/'!,.) = -e^^'^^^/'^X^^a. 

If we replace e'^'^™^/'V by f, then the solvability result Proposition 12. 11 shows that 
we can find a solution f, so that a solution r to (|2.13p is given by r = e~" '^™?/^f. 

To get a norm estimate for r, notice that for the right hand side of (|2.13p we 
have 

as /i ^ 0. The solvability result E?T] gives then that 

lkli/L(n)=0(/^'/'), 

as h 0, which implies that Hj^Hz/i ^(n) = 0{h^/^), as ft. — ?> 0. 

Thus we have obtained the following existence result for complex geometric optics 
solutions. 

Proposition 2.4. Let A &W'^^°'={VL,R^) and q e L°"{n,C). Then for h > small 
enough, there exist solutions u G H^{il), of the equation 

La.qU = in fi, 

that are of the form 

u{x,C:h)=e-</\a{x,C,h)+r{x,C,h)), 

where C G C^, is of the form given by (j2.6p . a G C°°(ri) solves the equation (j2.8p . 
and where a and r satisfy the estimates 

II alUoo(o) < C„ft-l"l/3 and Mn^^^in) = 0{h^'^). 

□ 

Remark 2.5. /n i/ie sequel, we need complex geometric optics solutions belonging 
to H'^ifl). To obtain such solutions, let fl' ft be a bounded domain with smooth 
boundary, and let us extend A G R^) andq€ L°°{n) to W^^°°{n',R^) and 

L°°{fl')- functions, respectively. By elliptic regularity, the complex geometric optics 
solutions, constructed on Q,' , according to Provosition 12.41 belong to H^(fl). 

Remark 2.6. Recall that <^ = ^N-'^{-i{a+il3)-{A'ioT~'^))oT. LemmaWBimvlies 
that : Co(il) — > C{Vl) is continuous. The estimates (j2.5p show that A^ ^ A 
uniformly onil. It follows that, if we define := ■^N^^{—i{a + if3) ■ {AoT^^))oT , 
then solves the equation 

(2.14) Co ■ V$" = -<o -A in fl, 
and satisfies 

||$(x,Co;ft'/')-<i>°||L-(o) ->0, ft->0. 



7 



Remark 2.7. We shall later use a slightly more general form for the amplitude a 
in the complex geometric optics solutions. Namely we suppose that a ~ ge*, where 
g £ C°°{^1), is such that 

(2.15) Co • Vg = 0. 

This means that g is holomorphic in a plane spanned by a and /3. Notice also that 
by picking a = ge*, we get by (|2.8p that 

Co-gV$ = -<o-.9^*, 

in place of (|2.9|1 . But the $ solving (j2.9p also solves the above equation. Hence we 
can use the same argument to obtain the $ for the above equation, as earlier. 
We thus obtain CGO solutions of the form 

u = e^</\ge'' + rg), 

where <f> solves (j2.9D . 

Notice also that setting a — ge^ does not affect the norm estimates on a in 
Proposition \2.4[ since g does not depend on h. 

3. An integral identity 

One central step in the ideas that are used in proving uniqueness results for 
inverse boundary value problems, is to derive an integral equation that expresses 
LP' orthogonality between the product of two solutions ui and 1*2, and the difference 
of two potentials qi and (72, see [23]. One shows that 

j [qi - q2)uiU2 = 0, 

provided that the DN-maps for qi and q2 are equal. 

A similar thing will be done in this subsection, for the magnetic case. The 
integral equation, is however more involved in the case of a magnetic potential and 
will not by itself be interpreted as an orthogonality relation. We will be considering 
the integral equation in conjunction with solutions that depend on a small positive 
parameter h. In the later sections we will see that in the limit /i — >■ 0, we obtain a 
criterion for the curl being zero. 

It will be convenient to set 

I dR^nN, 

Recall that we assume that 

We can thus choose Tj, such that 

fj C Tj, fj CC dRl \N, j = 1, 2. 
Then it follows from (jl.6p that 

(3.1) A^,,,,(/)|p^ =A^,,,,(/)|p^, 

for any / £ _ff'^/^(9R^), supp(/) C r2. In order to prove Theorem 11.21 we shall 
only use the data (|3.ip . which turns out to be enough to determine the magnetic 
field and the electric potential. 

We now begin deriving the integral identity. We assume that Aj,qj and Tj are 
as in Theorem [12] so that (|3.1|) also applies. 

Let ui G i7]Qp(]R'^) be the radiating solution to 

{LA,.,g,-k^)ui^OmRl, 
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with / e H'^^'^idRl), supp(/) C f 2. Let w e H^^^{W_) be the radiatmg solution to 

(La2,92 - k^)v = in RI, 

Define w :— v — ui. Then 

(La2,92 - k'^)w = 2i{A2 - Ai) ■ Vui + iV ■ {A2 - Ai)ui 

(3.2) + {Af ^ Al)ui + {qi - q2)ui. 
It follows from p.ip that 

(3.3) (5,1 + iAi ■ n)ui|p^ = (a„ + iA2 ■ n)v\f.^. 

By Lemma [1.11 we may and shall assmne that Ai ■ n = A2 ■ n = on dM.^_, so 
that dnW ~ on Fi. We also conclude from p.2p that w satisfies the Helmholtz 
equation 

{-A - k^)w = in Rl\N. 
As iflp^ = dn w|pj = 0, by unique continuation, we get that w = in R'^ \ N (see 
Theorem 16.31 and Corollarv l6.4l in the appendix). Since w E H'^^^{M.^_), we have 

w = d„w = on dNnW[. 

Let U2 € H'^{N) be a solution to (La2,w ^ k'^)'U2 = in N. Then by Green's 
formula, we get 

((-^^^2,92 - k^)w,U2)L'^(N) = (w, (iA2,W - ^^)"2)L2(Ar) 

- ((9,1 + 1^2 • n)w, U2)L^(dN) 

+ (w, (i9„ + iA2 ■ n)u2)L^(dN) 

= -(a,jW,U2)L2(,). 

Assuming that 

M2 = on I, 

we conclude that 

((^A2,g2 - fc^)w,U2)L2(Ar) = 0. 

Using equation (j3.2p we may write this as follows, 

(2i{A2 - Ai)- (Vui)iIJ + iV • (A2 - Ai)uiu^) dx 

+ {A\~ Al + qi- q2)uiu^dx = 0. 

JN 

Using again the fact that {A2 — Ai) ■ n ~ Q on d N and an integration by parts, we 
get 

if V • {A2 — Ai)uiU2dx = —i {A2 — Ai) ■ (VuiU2 + ui\7u2)dx. 
Jn Jn 

Thus, we obtain that 

/ i{A2 — Ai) ■ (VuiU2 ~ W1VU2) dx 
Jn 

+ / {Aj-Al + qi- q2)uiu^dx = 0, 

where m € Wi(M.l) and U2 e WHN). Here 

Wi{R^_) := {u e Hl^iW) : {LA,,q, - k^)u = in M.^_, 
supp(ui|aj{3_) C r2,u radiating}. 
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and 

W;{N) -.^ {u£ H^{N) : {LA^,q^ -k^)u^O in N,u\i = 0}. 

We shall next extend the integral identity p.4p to a richer class of solutions to 
the magnetic Schrodinger equation. To that end, let us introduce the following 
space of solutions, 

Wi{N) ~ {u e H^{N) : {LA,.q, - k^)u = in N,u\i = 0}. 

The following Rungc type approximation result is similar to those found in [S] , |16) 
and [TT]. 

Lemma 3.1. The space V := Wi{R^_)\n is dense in Wi{N) in the L'^{N) -topology. 

Proof. Suppose that V is not dense in Wi{N). First notice that span{V) = 1^ so 
that y is a linear subspace of L?{N). Since V is not dense in Wi{N), we have a 
vector uq S Wi{N) such that ^ V . We can decompose uq as uq = a + fe, where 

a £ V , b £ and 6 7^ 0. Let T be the linear functional on L'^{N), defined by 
T{x) := j>Toiy±(x)/\\b\\L2, where proj.^^ is the orthogonal projection to V . Now 
clearly \\T{uo)\\l^- = 1 and T\v = 0. 

By the Riesz representation theorem, there is gr G L^{N) that corresponds to 
T. Extend gr by zero to the complement of N in M"^. Let U G 7J[^^(R'^) be the 
incoming solution to 

{LA,,gT - k^W = gr in 
U\dMl = 0. 

The existence of such a solution is proved in [19] . 

Now let u e W^i(R^). Then because T\v = and supp(gT) C N, we get by the 
Green's formula of Lemma 16.21 that 

= (u,5t)l2(r3j = (m, {LA^^qT - A:^)C/)i2(H3j 

= ((-^^1,91 - [/)l2(k3j 

- {u, {dn + iAi ■ n)C/)i2(aH3_) 
+ ((9„ + iAi ■ n)u, f7)L2(aM3j 
= -(u,a„?7)^2(f,). 

Since the boundary condition u|p^ can be chosen arbitrarily from C^{T2), we get 
that 9„[/|p^ — 0. Since U\p^ = 0, we apply the unique continuation principle to 

conclude that U\^3 = 0. As [/ e iJjQ^(R^_), we have 

U\dNnS.^_ = 9n C^|aArnR3_ = 0. 

Now applying Green's formula and doing the same computation as above for uq 
and TV instead of u yields 

- {uq, {dn + iAi ■ n)U)L2(QN) 
+ {{dn + iAi ■ n)uQ, U)l2(on) 

= -{uQ,dnU)L2(l) = 0. 

Here we have used that uq\i — 0. It follows that T{uo) = 0. This contradiction 
completes the proof. □ 
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Since {A2 — Ai) ■ n ^ on d N, we can rewrite p.4p in the following form, 

- / uiiV ■ {{A2 ~ Ai)u^) dx- i{A2 - Ai) ■ (wiVmJ) dx 
Jn Jn 

+ / (Aj-Aj+qi- q2)uiu^dx = 0. 
J N 

Hence, an application of Lemma l3. II implies that the integral identity p.4p is valid 
for any ui € Wi{N) and any U2 € Wi{N). 

We summarize the discussion in this subsection in the following result. 

Proposition 3.2. Assume that Aj,qj andVj, j = 1,2 are as in Theorem M.'A and 
that the DN-maps satisfy 

(3.5) A^,,,,(/)|r, -A^,,,,(/)|r,, 

for any f g Hlllp{dM?_), supp(/) C Tj. Then 

I i{A2 — ^1) • (VwiU2 — uiWui) dx 

(3.6) 

+ / {A\- Al + qi-q2)uiu^dx = 0, 



IN 

for any ui G Wi{N) and any U2 G W2{N). 

□ 



Remark 3.3. Notice that the proof of Provosition \ 3.2\ only uses the assumption 
p.ip , which follows from p.Sp . Proposition Ig.gl holds therefore also under the 
weaker assumption p.ip . 



4. Recovering the magnetic field 

The aim of this section is to prove the first part of Theorem ll.2l by showing that 
the curl of the magnetic potential is determined by the DN-map. We choose an 
open ball B centered on dM.^_ with N CC B. And use the notations 

B_:=RlnB, B+:=M^nS, lB:=dRlnB. 

The first step in the argument will be to construct complex geometric optics so- 
lutions ui G Wi{N) and U2 G VF2*(iV) and then to examine the limit of p.6p as 

For Ml G Wi{N) and U2 G W2*(A^), we have that u^\i = 0, j = 1,2. To obtain 
solutions that satisfy this condition, we will first choose solutions defined on the 
larger set B and then use a reflection argument. 

The parameters C for the complex geometric optics solutions will be picked as 
follows. We will assume that 

(4.1) ^,71:72 e M^, I71I = I72I = 1 and that {71,72,^} is orthogonal. 

Similarly to [22], we set 



(4.2) Ci = ^+7i+*yi-/i^^72, 



C2 = ^ - 71 + «Y 1 - 72, 

so that Q ■ Q = 0, j = 1, 2, and (Ci + ^2)/^ = i^- Here ft. > is a small semiclassical 
parameter. Notice also that (j, j = 1,2 satisfy the conditions on ( in (|2.6p . when 
we take a = ±71 and /? = 72. 



11 



We need to extend the potentials Aj and gj, j = 1, 2, to _B+. For the component 
functions Aj^i, Aj^2, and qj, we do an even extension, and for ^j,3, we do an odd 
extension, i.e., for j = 1,2 we set, 

A, = h^'^'^^^ k^l,2, 



qj{x) 



Aj,3{x), X3 < 0, 

^-Aj^x), X3 > 0, 

qjix), X3 < 0, 
qjix), X3 > 0, 



where x := {xi,X2, —x^). By our assumptions, AjAx3=o — 0, from which it follows 
that Aj e and q^ € L°°{B), j = 1,2. 

We can now by Proposition 12.41 and Remark 12.51 pick complex geometric optics 
solutions ui in H^{B), 

uiix,Ci;h) = e^<i/''(e*i(^'''i+'''^''') +ri(2;,Ci;/^)) 

of the equation (L^^ - k'^)ui = in B, where $i £ C°°(B). By Remark 
$1 — >■ $5 in the L°°-norm as /i 0, where $° solves the equation 

(4.3) (7i+i72)-V$? = -i(7i+i72)-ii in S. 
To obtain a function that is zero on the plane 3:3 = 0, we set 

(4.4) ui{x) ui{x) ~ ui{x), xeB^UIb- 

Then it is easy to check that the restriction ui\n € Wi{N). 

We can similarly pick by Proposition 12.41 and Remark 12.51 complex geometric 
optics solutions 112 in H^{B), 

Mx, C2; h) = e"<^/''(e*=("'-^i+'^^''') + 7-2(2;, C2;h)) 

of the equation (L^^ - k'^)u2 = in B, where $2 £ C°°(B). By Remark 
$2 ^2 in the L°°-norm as /i 0, where $2 solves the equation 

(4.5) (-71 + ^72) ■ = -i(-7i + 172) ■ A2 in B. 
To obtain a function that is zero on the plane X3 = 0, we set 

(4.6) U2{x) := U2{x) — U2{x) , xQB^UIb- 

Then it is easy to check that the restriction U2\n ^ 1^2* (-^)- 

The next step is to substitute the complex geometric optics solutions ui and U2, 
given by ()4.4p and (j4.6|) . respectively, into the integral identity (j3.6|) . This will be 
done in the Lemma bellow. We will use the abbreviations ^1(2;) := e*^^^^ + ri{x) 
and P2(a;) := e*"^'') + r2{x), so that 

Mi(x) = e^<i/''Pi(x) - e*<i/''Pi(i), 

^2(2;) = e'=<^'^P2{x) - e^<^'^P2{i). 

For future references, it will be convenient to compute the product of the phases 
that occur in the terms uiu2, VwiU2 and uiVu2 

(4J) ^i-Ci/h^x^/h _ ^ix-(,^i{0,0,-2x3)<i/h _ ^ix-(,--2-ii^3X3/h ^ 

gX-Cl/h^x-Q-z/h _ gi£-4gj(0,0,22;3)-Ci/h _ ^ix-(,+ +2^i^3X3 / h 
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where 7^- = (7j,i> 7j,2, 7j,3), j = 1,2 and 



?± = (^1'^2,±-^y1 ^72,3 

We restrict the choices of 71 and 72, by assuming that 
(4.8) 71,3 = and 72,3 + 0. 

We need these conditions for the proof of the next Lemma. The first condition 
makes the above phases purely imaginary, which avoids exponential growth of the 
terms, as /i — > . The second condition implies that |^±|— T-ooasft.— ?>0. This will 
be needed since we will use the Riemann-Lebesgue Lemma to eliminate unwanted 
imaginary exponentials. 

Finally it will also be convenient to explicitly state the following norm estimates, 
which follow from Proposition [ 



||e*^|U. =0(1), |!Ve*^|U~=0(/.-i/3), 
\\r.^\^.=0{h}l^), \\^T,\\^.^0{h-^l\ j = l,2, 

as h 0. 

Lemma 4.1. We have 

(4.10) (71 + 172) ■ f {A2- ii)e"-«e*?+*^dx = 0, 



where 71 , 72 and ^ satisfy (|4.1|) and (j4.8l) . 

Proof. Wc will prove the statement by multiplying the integral equation (j3.6p of 
Proposition 13 . 21 bv h, when ui and U2 arc given by (j4.4p and (j4.6p . and then take 
the limit as ft, — !• 0. 

To begin with notice that we may integrate over B_ in p.6p . since 

supp(Aj), supp((3'j) C N d 

and Uj arc defined in i?, when j — 1, 2. We first show that for the second term in 
p.6p we have 



(4.11) h {Aj - Al + qi - q2)uiU2dx ^ 0, 

Jb 



as /i ^ 0. Using the phase computations (|4.7p we get that 

u,u2 = e^«-"Pi(x)P2(x) - e"-«+Pi(x)P2(i) 
- e"^-«-Pi(i)P2(x) + e'«-^Pi(x)P2(i). 

This is multiplied by an L°° function in (|4.1ip . Since we restricted the choice of 
71 to make the exponents purely imaginary, we see easily using the estimates (j4.9p 
that (|4TT|) holds. 

Equation (|3.6p multiplied by h, is thus reduced in the limit to 



(4.12) lini ^/i J i{A2 - Ai) ■ VuiU2dx j i{A2 - Ai) ■ ui'\/u2dx^ = 0. 
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We will proceed by examining the first term. Using (j4.7p we write VuiM2 as 

+ e'^'-^VPiixyP^ - e"-«+VPi(x)^^(I) 

_ ^(e" «-Pi(i)/M^- e" «Pi(i)7M5)) 
- e'^-«- V,Pi + e"-« V,Pi , 

where Q := (Cj,i, Ci,2j ~Cj,'i)i j = I7 2. The terms of the product that do not contain 
the factor 1/h, result in integrals similar to the one in (j4.1ip . One sees similarly 
using estimates (|4.9p that they are zero in the limit of ()4.12[) . The first term inside 
the limit in (j4.12p is therefore reduced to 

lim / i{A2 - Ai) • (Cie" «Pi(x)/M^ - Cie" «-Pi(5)/M^ 

-Cie"-«^ Pi {x)lW) + Cie"-«Pi {i)lW))dx. 

Now we use the Riemann-Lebesgue Lemma to conclude that the terms with ex- 
ponents containing ^-|_ and ^_ arc zero in the limit. To see this, notice that by 
Remark 1 2. 6 1 we see that ||^'i||L°°(B_) < for some C > 0, when h is small enough. 

Estimates (g^l) show that \\r^\\L^(B ) = 0{h^/'^). Hence ||PjPj||li(b_) < C, for 
some C > when h is small enough. Finally we have ^± ^ 00, as /i — > 0, because 
of the restrictions (|4.8p . 

The first term in (j4.12p is therefore 

lim / i(A2 - Ai) • (Cie" «Pi(x)/M^ + Cie" «Pi(i)i^)da; 

as ft, — > 0. The terms containing ri in the products of Pi and P2 are, because of 
(j4.9p . zero in the limit. The above limit is thus equal to 

lim / i{A2 - Ai) ■ (Cie"-«e*i(")+*^ + Cie"-«e*i(^)+*^)dx. 

Finally we split the integral and do a change of variable in the second term and 
arrive at the expression 

(4.13) lim / i(i2 -ii) •Cie""^e*i(")+*^da;, 

for the first term of (|4.12|) . 

Returning to the second term in (|4.12p . containing uiV?l2. This is of the same 
form as the first one. By doing the above derivation by simply exchanging the roles 
of ui and M2, we similarly see that the second term becomes 

(4.14) lim - / i(i2-ii)-^e"-«e*i(^)+*^da;. 

Now Ci — > (71 + ij2) and ^2 — (71 + «72), as h ^ 0. Thus by using (|4.13p with 
(j4.14p . we can rewrite (|4.12p as 

lim / iiA2 - ii) • (Cie"-«e*i(")+*^ _ <2e"-^e*i(")+*^)(ia; 
= 2 / t{A2 - ii) • (71 + i72)e"-«e*i'(")+*^da; = 0. 

JB 

□ 
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The next Proposition shows that (|4.10p holds even when the exponential function 
depending on z = 1,2 is removed. The argument follows [5] closely. We will 
give details for the convenience of the reader. 

Proposition 4.2. The equality (j4.10p implies that 

(4.15) (71 + *72) ■ f {A2- ii)e"-«(ix = 0, 

for 7i , 72 and ^ which satisfy (|4.ip and (|4.8p . 
Proo/. By and we have that 

(4.16) (7i+i72)-V($; + $f) = -i(7i + i72)-(ii-^) in B. 

Remark 12 . 71 furthermore implies that the amplitude e*^ in the definition of wi can 
be replaced by ge"^^, if g G C°°{B) is a solution of 

(4.17) (71 +«72) ■ V.g = in B. 

Let *(a;) := $?(a;) + $2(0;). Then instead of (|4.10p we can write, 

(7i+«72)- / (i2-ii)5e"-«e*(")dx = 0. 

We conclude from (|4.16|) that 

(71 + 172) • (^2 - ^i)5e* = -«(7i + «72) • (fiVe*), 
and therefore, we get 

(4.18) / 5e"«(7i+i72)-Ve*da; = 0, 



for all g satisfying ()4.17p . 

We pick a 73, with I73I = 1, so that we obtain an orthonormal basis {71, 72, 73}- 
Let T be the coordinate transform into this basis, i.e. y = Tx = (a; ■ 71 , a; • 72, a; • 73). 

Set z ^ yi + iy2, so that dz = [dy^ +i dy^)/2 and 

(71+472) •V = 29j. 
Rewriting (j4.18p using this and a change of variable given by T we have 



JTB 



ITB 

for all g satisfying (|4.17p . 

Notice that y = y^S,^, since ^ is in the y-coordinates of the form (0, 0, ^3). The 
above integral is therefore a Fourier transform w.r.t. ^3. Let g € C°°{TB) satisfy 
dz 9 = and be independent of t/3. Then taking the inverse Fourier transform we 
write 



= / gd-e'^dyidy2 
d-{ge^)dyidy2, 



where Ty^ := TB n IIj,, and Wy^ ~ {(2/172/2,2/3) : (2/1,2/2) G R^}- Notice that the 
boundary of Ty.^ is smooth. Multiplying the above by 2i and using Stokes' theorem 



15 



we get that 

= 2i d-{ge^)dyidy2 

JTy^ 

= / V X (.ge*,ige*,0) • ndyidy2 
= / (5e*,*ge*,0)-dZ 

JdTy^ 

(4.19) = / .ge*dz, 



for all holomorphic functions g G C°°{Ty^). 

Next we shall show that (j4.19p implies that there exists a nowhere vanishing 
holomorphic function F € C{Ty^) such that 

(4.20) F\9Ty,=e''\aTy,. 
To this end, we define F to be 

1 f e*(«) 

F(.) . - y^^^^^ ^dC, .eC\dTy^. 

The function F is holomorphic away from dTy^. As e* is Lipschitz, we know 
because of the Plemelj-Sokhotski-Privalov formula (see e.g. [10]), that 

(4.21) lim F{z)- lim i^(z) = e*(^°\ zoGdTy^. 

z-s-zo,zGT„3 z->zo,z^T„3 

Now the function C (C ~ -2^)"^ is holomorphic on Ty^ when z ^ Ty^. By choosing 
(;(^) = (( - z)"^ in (|4.19p . we get therefore that F{z) — 0, when z ^ Ty^. Hence, 
the second limit in (j4.2ip vanishes, and therefore, F is holomorphic function on 
Ty3, such that holds. 

Next we show that F is non- vanishing in Ty^. When doing so, let dTy^ be 
parametrized by z = j(t), and N be the number of zeros of F in Ty^. Then by the 
argument principle, we get 



27rt F{z) 2-Ki J Pay ( 2Tri J^^o-, ( 

To see that the last integral is zero, notice that this the winding number of the 
path 6*°'''. And that e*^'''^*)^ is homotopic to the constant contour {1}, with the 
homotopy given by e**'^''''*-'-' , s G [0, 1]. 

Next, since _F is a non- vanishing holomorphic function on Ty^ and Ty^ is simply 
connected, it admits a holomorphic logarithm. Hence, (|4.20p implies that 

(log^^)|aT„3 = 1'|aT„3. 
Because log_F = is continuous on dTy^. wc have by the Cauchy theorem, 



g'i'dz ^ / ghgFdz^O, 

dTy, JdTy, 

where g £ C°°{Ty^) is an arbitrary function such that d^g ~ 0. Using Stokes' 
formula as in (j4.19[) allows us to write this as 



gds'^dyidy2 = 0. 
Taking the Fourier transform with respect to y^ , we get 



JT(B) 
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for all ^ = (0,0,^3), ^3 G M. Hence, returning back to the x variable, we obtain 
that 

{11+112)- I e"-«5(x)V*(a;)dx- 0, 

J B 

where g e C^(B) is such that (71 + ^72) • V.g = in B. 
Using (|4.16p . we finally get 

(4.22) (71 + 172) • / (is - Ai)g{x)e'-<dx = 0. 
Setting 5 = 1, we obtain (|4.15p . 

□ 

By replacing the vector 72 by —72 in (j4.15p . we see that 

(4.23) (7i-«72)- / (^2 -ii)e"-«(ix = 0. 

Jb 

Hence, and imply that 

(4.24) 7 • / (^2 - ii)e" «d.T = 0, 

JB 

for aU 7 G span{7i,72} and aU ^ G such that and hold. 

In the proof of the next Proposition we see that (|4.24D is actually a condition 
for having V x (Ai — A2) =0. This is therefore the last step in proving that the 
DN-map determines the curl of the magnetic potential. 

Proposition 4.3. It follows from (j4.24p that 

(4.25) V X ii = V X ^2 in B. 

Proof. Assume that ^ G M'^ is not on the line L := (0, 0, i), t G M. Then the vectors 
71 and 72 given by 

71 := (-6,6,0), 7i:=7i/|7i|, 

(4.26) 72:=Cx7i, 72:= 72/|72|, 

where ^ x 71 stands for the cross product, satisfy (|4.8p and ()4.ip . Thus, for any 
vector i G \ L, (|4:24)) says that 

(4.27) 7-^^(0-0, v{0:=A^xiO~^x{0, 

for all 7 G span{7i, 72}. Here x is the characteristic function of the set B. For any 
vector ^ G M^, wc have the following decomposition, 

where Re 11^(6, Iin?;^(C) are multiples of ^, and Rei'_L(6, Im?^_L(C) are orthogonal 
to ^. Now we have Re (6 , Im w_l (6 G span{7i,72}, and therefore, it follows from 
(irt?]) that v±{^) = 0, for aU ^ G \ L. 
Hence, v{^) = so that that 

e X v{0 = 0, 

for all ^ G R'^ \i, and thus, everywhere, by the analyticity of the Fourier transform. 
Taking the inverse Fourier transform, we obtain (|4?25l) . □ 
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5. Determining the electric potential 

In order to complete the proof of Theorem 1 1.21 we need to show that the electric 
potential is also determined by the DN-map. Again we assume that Aj , qj and Tj , 
j = 1,2 are as in Theorem ll.2l and that the DN-maps satisfy (|1.6p . and hence p.ip . 

Since B is simply connected, it follows from the Helmholtz decomposition of 
Ai ~ A2 and (|4.25p that there exists ip G C^'^(B) with i' = near dB such that 

Ai = A2 + Vi/j in B. 

We extend ip to a function of class C^'^ on all of M"^ such that tp = on M.^ \ B. 
Then 

ii=A2+V-0 in K'\ 

Since Tj C R^\N, j ^ 1,2, and \ Iv is connected, we can assume that "0 = 

on M.'i \ N and hence, we have that -0 = on Fi U It follows then from Lemma 
ll.ll part (i) and p.l[) that for all / with supp(/) C 

Therefore, we can apply Proposition 13.21 with Ai = A2 and get 

(5.1) / (gi - q2)uiu^dx = 0, 

Jn 

for aU ui £ Wi{N) and U2 E W^{N). 

Choosing in (|5.ip ui and U2 as the complex geometric optics solutions, given by 
(|4.4p and ()4.6p . passing to _B_, and letting h ^ 0, we have 

(5.2) / {qi - q2)e"-«e*?(")+*^dx = 0. 
Jb 

By RemarklOe^i in the definition (031) of ui can be replaced by .ge*i if g G C"^ (i?) 
is a solution of 

(71 + 172) • V5 = in 
Then (j5.2l) can be replaced by 

/ (gi - q2)3(x)e"-«e*°(")+*^dx = 0. 
Jb 

Now (|4.16p has the form, 

(7i+i72)-V($? + $0) = in B, 
since we have that Ai — A2. Thus, we can take g = e^'*i+*2) and obtain that 

(5.3) / (gi - g2)e" «dx = 0, 

Jb 

for all ^ e such that there exist 71,72 € K^, satisfying (|4.ip and (j4.8p . Since 
for any ^ e R'^ not of the form ^ = (0,0,^3), the vectors, given by (j4.26p . satisfy 
(|iTl) and we conclude that ([O)) holds for all ^ G R^ except those of the form 

5 = (0, 0, ^3), and therefore, by analyticity of the Fourier transform, for all ^ G R^. 
Hence, qi = q2 in S_. This completes the proof of Theorem 11.21 
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6. Appendix 

6.1. Magnetic Green's formulas. Let us first recall, following [5], the standard 
Green formula applied to the magnetic Schrodinger operator. 

Lemma 6.1. Suppose that C R'^ is open and bounded, with piecewise bound- 
ary. Let A e M^i'°°(rj,K3) and q e L°°{n). Then we have, 

= (ti, {dn + iA ■ n)v)L'2(Q^) - {{dn + iA ■ n)u,v)L2(Q^), 
for all u,v e H^{Vl), with Am, Az; G L'^{VI), where n is the exterior unit normal to 

dn. 

We shall also need a version of the above result where fl is replaced by M?. . We 
shall then need to put some restrictions on v and u, because R'l is unbounded. To 
this end wc assume that u and v arc solutions to the Hclmholtz equation outside 
some compact set, that obey some form of radiation condition. To be precise, let 
A e W^i^p{Rl,R^), q e LS™p(]R^-), and let u e Hl^{W_) be such that 

{LA^q~k^)u^O in M.^, 

supp(u|g][{3 ) is compact, and u is outgoing. Assume also that v £ -f^iocC^-) satisfies 

supp(w|gjj3 ) is compact, and v is incoming. 
Lemma 6.2. With u and v as above, we have 

{{La.q - k'^)u, w)l2(r3j - {u, {LA,q - fc^)w)L2(R3j 

^ {u, {dn+iA ■ n)v)L2(aR3j - {{dn + iA ■ n)u, u)i2(gR3j. 

Proof. Let Bji := B{xo,R) be an open ball in M'^ of radius R, and choose R > 
large enough so that 

supp(A),supp(g) C Br. 
Set Q = R^_ n Br. By Lemma [Ql we know that 

{{LA,q - k^)u, v)L2(f2) - {u, {LA,q ~ fc^)w)L2(n) 

= {u, {dn + iA ■ n)v)L2(Q^) - {{dn + iA ■ n)u,v)L2(Qa)- 
Thus, to obtain (|6.ip we need to show that 

(6.2) / {udnV — dnUv)dSR, — > 0, i? — oo. 

Let us rewrite the left hand side of the above as follows, 

/ {dnV — ikv)udSR — / (9„u — iku)vdSR. 

We show that first term goes to zero as i? cxd. The second term can be handled 
in the same way. Applying Cauchy-Sehwarz gives 

2 



{dnV — ikv)udSR, 



< I \dnV~lkv\^dSR f H'^dSR. 



Here the first integral goes to zero, since dnV — ikv = dnV + ikv and \dnV + ikv\'^ 
is o(l/r^) as r = |a;| — )■ cx), since v is incoming. We conclude the proof by showing 
that the second integral is bounded as i? oo. 



19 



To this end we let R2 > Ri, where Ri is such that 

supp(A),supp((7) C B(.To,i?i), 

and set Bj := B{xo,Rj), j = 1,2 and U {B2 \ Bi) n M?_. We multiply the 
Sommerfeld condition (|1.4p by its complex conjugate and integrate over d B2 flM'^, 
which gives 

(fc2|u|2 + \^r^u\'^ +2klm.{udnU))dS 

d B2 nKl 

(6.3) = / \dnU~iku\^dS ^0, 

J d B2r\V?_ 

'AS R2 ^ 00 and where n is the outer unit normal vector to B2. 
By Green's formula we have on the other hand that 

(6.4) / udnU — udnU = / uAu — Auu = 0. 
Jdu Ju 

We may assume that wlgj/^g^s = 0, since supp(w|9]R^) is compact. We can thus 

write (16.41) as 



(6.5) / Im(ti9„u) = / Im(7i9„w). 

But this implies that the / |wp and /|9„up terms stay bounded in the limit 
(lOl) . □ 

6.2. The unique continuation principle. In this work we make heavy use of 
the so called unique continuation principle. The unique continuation principle can 
be seen as an extension of the familiar property that an analytic function that is 
zero on some open set is identically zero. 

Let il C M" be an open connected set, and let 

Pu = o^iji^) didjU + bi{x) diU + c{x)u. 

i,j=l i 

Here Oij G C^{fl) are real- valued, Oij = aji, and there is C > so that 

n 

Furthermore, bi € L°°{^l,C) and c G L°°(il,C). We have the following result, see 
[g and [H]. 

Theorem 6.3. Let u G H^^^{Q) be such that Pu = m 17. Let ui C fl be open 
non-empty. If u = on lo, then u vanishes identically in fl. 

Corollary 6.4. Assume that dO, is of class . Let T CZ dfl be open non-empty. 
Let u G H'^{il) be such that Pu = in fl. Assume that 

u — BuU = on r. 

Here B^u is the conormal derivative of u, given by 

n 

i.j=i 

Then u vanishes identically in ^l. 
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